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According to recent reports there is an excess correlation and an apparent regularity in the galaxy 
one-dimensional polar distribution with a characteristic scale of 128 h" 1 Mpc. This aparent spatial 
periodicity can be naturally explained by a time oscillation of the gravitational constant G. On 
the other hand, periodic growth features of bivalve and coral fossiles appear to show a periodic 
component in the time dependence of the number of days per year. In this letter we show that a 
time oscillating gravitational constant with similar period and amplitude can explain such a feature. 
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Our theoretical description of nature presents a set of parameters, the so called fundamental constants, that have 
to be determined from experience. It is generally thougth that this is so because we lack a unified theory of all 
interactions. The time variation of any fundamental constant might be one of the few low energy phenomena which 
could manifest this "new physics" . The time evolution of these parameters is supposed to be governed by dynamics 
of cosmological origin, so its variation rate is likely to be of the same order (or less) than the Hubble rate. We should 
seek then for low energy and long term (millions of years) phenomena pH§ or else for high precision measurements 
(ten or more places in the decimal expansion) of some "constant" parameter, separated by a few years |§J^]. Actually, 
superstring theories |j) and Kaluza-Klein theories have cosmological solutions in which the low-energy fundamental 
constants vary with time [8|-[l2| . 

Broadhurst, Ellis, Koo and Szalay [ p"3[ combined data from four distinct deep pencil-beam surveys at the north and 
south Galactic poles to produce a well sampled distribution of galaxies by redshift on a linear scale extending to 2000 
h^ 1 Mpc. They reported a periodicity in the galaxy distribution of 128 h~ x Mpc. Soon afterwards Morikawa |Tj ] 
noted that this aparent spatial periodicity could be naturally explained by a time oscillation of the Hubble parameter. 
In his model the oscillation was produced by a massive scalar field non-minimally coupled to gravity, inducing also 
the time oscillation of the gravitational constant G. Hill, Steinhardt and Turner |0| proposed different scenarios, 
including a time oscillating Hubble parameter, but also the oscillation of atomic lines as an alternative explanation 
of the red-shift galaxy distribution. An oscillation in the Rydberg constant due to the variation of the fine structure 
constant a or the electron mass m e requires a modification of the standard model in which a or m e become dependent 
on a scalar field. Thus an oscillating a introduces a Yukawa potential between samples with non-zero electrostatic 
enegy contribution, while an oscillating m e induces a similar interaction now proportional to the lepton numbers of 
the sample s [fl6| . In this context both possibilities were shown to be ruled out jlfj] by the experiment of Braginsky 
and Panov |17| , leaving the Hubble oscillating scenario as a solely candidate to analyze. 

Even if the periodicity reported by Broadhurst et.al. |l3[ ] were an artifact of small angle sampling of galaxy distri- 
bution, the problem of an oscillating gravitational constant is interesting enough to deserve a different observational 
analysis. There are indeed superstring models that predict oscillating fundamental constants. These are functions of 
massive dilaton fields associated with the size of extra dimensions @||- To account a period ct ~ \28h~ 1 Mpc the 
dilaton mass m must be O(10~ 31 eV). In Ref. the time dependence proposed for the gravitational constant is |l[| 

G = G + Gicos[2m(t-t )+ip], (1) 
and the equation for the cosmological scale factor is 

H> = -%H + (2) 

where H is the Hubble parameter. From Eq. |^ it is easy to see that the leading correction to the Hubble parameter 
is H ~ H[l — ^(G/GH)]. Then, since H = a/a and 1 + z oc a , one finds that 

dz G 

^ = l-^j = l+./lcos[2m(t-io)+# (3) 
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where we have defined the amplitude A — —mGi/HGo- Following Ref. |i5[ , if the Universe is spatialy homogeneous, 
with a uniform galaxy density per comoving volume no, the number of galaxies dN in a solid angle d£l with red shift 
between z and z + dz is modulated compared to the distribution in the absence of oscillations in the following way: 



dN 



2 dzdQ 



dN \ dz 
z^dz^d^l J dz 
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where to lowest order z 2 



We have then an apparent variation in the density of galaxies which is isotropic and 



has peaks lying on concentric spherical shells at periodically spaced radii. An amplitude A ~ 0.5 (or equivalently 
about 1% for G\/Gq), a phase ip ~ and a period cr ~ L = 128/i _1 Mpc approximately reproduces the result of Ref. 
|l3f . Such a period is similar to the time scales involved in the paleontological record. 

It is well known that several taxons record growth rhythms in their skeletons; i.e. periodic markings locked to the 
astronomical cycles of day, month and year Jl8|. From these growth rhythms, the number of days per year, days 
per month and months per year have been obtained as functions of geological time 18-2(3]. This parameters can be 
simply expressed in terms of the Earth's rotation fi, the mean motion of the Sun n Q and the mean motion of the 
moon n ) as follows: 



Nd/y = — 



TlQ 
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The effect of a varying G on the Earth-Moon system can be studied under the adiabatic hypothesis, as stated in 
references |2l],|22j , so the Keplerian equations of motion mantain their form, and G is replaced by the appropiate time 
function. Morover, in a reference system where G depends only on time but not on space, angular momentum is still 
conserved Ell] . Then we obtain 



n Q 



G' 
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where (3 ~ 1.83 depends on the mass and pressure distribution in the Earth's interior KM, 7 ~ 1.84 is related to the 
tidal couple of the Earth-Moon system |19|, h ) \ t is the tidal acceleration of the lunar longitude and —(2 — ($)% is the 
change of the moment of inertia due to the change of G 0]. According to Lambeck |2(]] the palaeontological data 
suggest that there may be significant systematic errors in the counts, and it is convenient to introduce parameters 
An0,An } andAfi accounting for such errors in the estimation of the planetary angular velocities. Then for a time 
oscillating G the final expressions for the observables are 
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N^ L _ l = _n l f+ An, + Ano (13) 
N", n. n. rip, 

d/m ) t ) U 

where i is minus the geological (not ephemeris) time plf and <fi = § — i ) - 

We used the data carefully filtered by Lambeck ]20| according to biological reliability criteria pj[ |. These data 
(37 points) are shown in Table |. Originally Lambeck adjusted the data assuming that only (constant) dissipation 
mechanisms are the responsible for f2 and n ) . The adjusted values of f2 and h ) were in good agreement with the modern 
astronomical values. The residuals showed no obvious systematic trends, indicating that the growth rings present a 
high degree of confidence. In our case we adjusted simultaneously the three curves with and without the oscillatory 
parameters (i.e. with and without an oscillatory G), using the Levenberg-Marquardt least-squares method. The x 2 



significances of the adjustments are shown in Table III. There is still a debate on the correct values for standard 
deviations of the palaeontological number counts, so the absolute value of the x 2 significance is not important. What 
it is really important is the change of the significance when we include the oscillatory parameters. We see that there 
is a conspicuous increase of significance when we include the oscillatory parameters. For this last adjustment the 



best fit values are shown in Table IV, together with independent 95% confidence limits. The upper confidence limit 
of Gi/Go is about 0.009, which is the desired amplitude accounting for the red-shift survey, i.e., we have marginal 
consistency with the oscillating G hypothesis. The best fit value of the period of oscillation remarkably coincides with 
the galaxy distribution period, although the confidence limits weaken its relevance. The values for (tidal) tl/ti . ,f2/f2 



and fiQ/nQ agree with other estimates [j20|. The phase ip is consistent with the zero value proposed in 15 1, which 
implies a zero value for the present rate of change of G. Our fit is then consistent with the current upper bounds 
on the time variation of G based upon the Viking radar-echo experiments [^| [note that only upper bounds based 
on present observations are valid if G oscillates, so the much more stringent upper bounds of |^l|,^2| do not apply 
because they are also based in long-term (several oscillation periods) phenomena]. In order to test the sensibility of 
the solution, we made the same adjustments using original data from several sources, as shown in Table [□], including 
bivalves, corals, cephalopods, brachiopods and estromatolites and totalizing 61 points. As shown in Table |l| in this 
case both significances are small, showing that the filtering of data as made by Lambeck introduces a bias towards 
the oscillatory hypothesis. 

Wc conclude that our results do not exclude an oscillating gravitational constant inducing a periodic galaxy distri- 
bution. Indeed the significance of the Lambeck data adjustment suggests thet there is an oscillatory component in the 
time evolution of planetary orbit ratios. However, there are several uncertainties in our model that forbid a definite 
conclusion. In the first place, paleontological growth rhythms are subject to large variations, and should be handled 
with care |23|-|25|l. Secondly, it should be noted that the changes in the resonance structure of the oceans due to 
continental drift provoke considerable variations of the Earth- Moon tidal torque within 100 million years time scales 
p6|. This fact could account as well for the time oscillation of the number of days per year. As a result we can only 
state that an upper bound for the oscillation amplitude of the gravitational constant is Gi/Gq < 0.01 (taken from the 
upper confidence limit of Table |i"v| ) . The Hubble oscillating hypothesis is strongly testable beacause it predicts the 
same periodic patterns in all directions. Also, as stated by Morikawa [|l4| both the density contrast and the distance 
to the nearest peak is not clear from the survey. Then it is worth making further deep-pencil surveys in the same and 
other directions in the sky. Summarizing, the coincidences which we observe in this work (coincidences may become 
consequences p7|) are significant enough so as to become a subject of further numerical and observational research. 



also at Departamento de Fisica, F.C.E.y N, Universidad Nacional de Mar del Plata, Funes 3350 - (7600) Mar del Plata, 
ARGENTINA 

[1] M. W. Mc.Elhinny, S. R. Taylor and D. J. Stevenson, Nature (London) 271, 316 (1978). 
[2] A. M. Wolfe, R. L. Brown and M. S. Roberts, Phys. Rev. Lett. 37, 179 (1976). 
[3] A. I. Shlyakhter, Atomki Report A/I (1983); Nature (London) 264 , 340 (1976). 
[4] T. Damour, G. W. Gibbons and J. H. Taylor, Phys. Rev. Lett. 61, 1151 (1988). 

[5] R. H. Hellings, P. J. Adams, J. D. Anderson, M. S. Keesey, E. L. Lau, E. M. Standish, V. M. Canuto and I. Goldman, 

Phys. Rev. Lett. 51, 1609 (1983). 
[6] J. H. Schwartz, Phys. Rep. 89, 223 (1982); M. B. Green J. H. Schwartz, Phys. Lett.l99B, 117 (1984). 
[7] T. Kaluza, Sitz. Preuss. Akad. Wiss. Phys. Math. Kl, 966 (1921); O. Z. Klein, Phys. 37, 895 (1926). 
[8] Y. Wu and Z. Wang, Phys. Rev. Lett. 57, 1978 (1986). 
[9] K. Maeda, Mod. Phys. Lett. A3, 243 (1988); 



3 



[10] S. M. Barr and P. K. Mohapatra, Phys. Rev. D38,3011 (1988). 

[11] A. Chodos and S. Detweiler, Phys. Rev. D21, 2167 (1980). M. Gleiser and J. G. Taylor, Phys. Rev. D31, 1904 (1985). 
[12] W. J. Marciano, Phys. Rev. Lett. 52, 489 (1984). 

[13] T. J. Broadhurst, R. S. Ellis, D. C. Koo and A. S. Szalay, Nature (London) 343, 726 (1990). 
[14] M.Morikawa, Astrophys. J. 369, 20 (1991). 

[15] C. T. Hill, P. J. Steinhardt and M. S. Turner, Fermi National Accelerator Laboratory prepint 90/129-T (1990). 
[16] D.Sudarsky, Fermi Institute preprint, 1992. 

[17] V.B.Braginskii and V.I.Panov, Zh. Eksp. Teor. Fiz. 61, 873 (1971). 

[18] C. T. Scrutton, in Tidal Friction and the Earth's Rotation (eds Brosche, P. Sundermann, J.) 154-196 (Springer, New York, 
1978). 

[19] K. Lambeck, The Earth's Variable Rotation (Cambridge University Press, New York, 1980). 

[20] K. Lambeck, in Tidal Friction and the Earth's Rotation, ed. by P. Brosche and J. Sundermann, (Springer- Verlag, Berlin 
1978). 

[21] P. D. Sisterna and H. Vucetich, Phys. Rev. D41, 1034 (1990). 
[22] P. D. Sisterna and H. Vucetich, Phys. Rev. D44, 3096 (1991). 
[23] R. A. Lutz and D. C. Rhoads, Science 198, 1222 (1977). 
[24] W. W. Hughes, Geophysical Surveys 7, 169 (1985). 
[25] G. D. Rosenberg, Geophysical Surveys 7, 185 (1985). 

[26] P. Brosche, in Sun and Planetary System, ed. by W. Fricke and G. Teleki (Reidel, D. Publishing Company). 

[27] J.D.Barrow and F.J.Tipler, The Anthropic Cosmological Principle (Oxford University Press, Oxford, 1988). 

[28] J.W.Wells, Nature (London) 197, 948 (1963). J.W.Wells, in Palaeogeophysics ed. by S.K.Runcorn (London, Academic 

Press, 1970). 
[29] C.T.Scrutton, Palaeontology 7, 552 (1965). 
[30] S.J.Mazzulo, Bull.Geol.Soc.Am. 82, 1085 (1971). 

[31] G.A.L. Johnson and J.R.Nudds, in Growth Rhythms and the History of the Earth's Rotation, ed. by G.D.Rosenberg and 

S.K.Runcorn (London, Wiley, 1975). 
[32] G.Panella, Astrophys.Space Sci. 16, 212 (1972). 

[33] W.B.N. Berry and R.M.Baker, in Growth Rhythms and the History of the Earth's Rotation, ed. by G.D.Rosenberg and 

S.K.Runcorn (London, Wiley, 1975). 
[34] W.B.N. Berry and R.M.Baker, Nature (London) 217, 938 (1968). 
[35] A.McGuhan, Abstr.Mtg.Geol.Soc.Am. 1967, 145 (1967). 

[36] S.M.Pompea, P.G.K.Kahn and R.B. Culver, Vistas in Astronomy 1979 23, 185 (1979). 



4 



TABLE I. Lambeck selection of counting and geological data. 



Geologic time [M.Y.] 


Counting 


Adopted S.D. 


Organisms 


Data type 


-300 


385 


4.2 


corals 


da 


-320 


398 


6 


corals 


da 


-370 


398 


1.7 


corals 


da 


-420 


400 


6 


corals 


da 


-440 


412 


6 


corals 


da 





360 


4 


corals 


da 


-370 


401 


6 


corals 


da 


-370 


30.66 


0.5 


corals 


dm 


-330 


30.2 


0.4 


corals 


dm 





359 


2.1 


bivalves 


da 


-70 


375 


3.5 


bivalves 


da 


-220 


372 


3.5 


bivalves 


da 


-290 


383 


3.5 


bivalves 


da 


-340 


398 


4.2 


bivalves 


da 


-360 


406 


6 


bivalves 


da 





29.2 


0.4 


bivalves 


dm 


-14 


29.4 


0.5 


bivalves 


dm 


-38 


29.8 


0.8 


bivalves 


dm 


-54 


29.6 


0.6 


bivalves 


dm 


-70 


29.9 


0.4 


bivalves 


dm 


-220 


29.7 


0.6 


bivalves 


dm 


-300 


30.2 


0.6 


bivalves 


dm 


-350 


30.4 


0.6 


bivalves 


dm 


-445 


30.3 


0.8 


bivalves 


dm 


-1 


14.75 


0.14 


bivalves 


dm/2 


-4 


14.83 


0.16 


bivalves 


dm/2 


-30 


14.82 


0.09 


bivalves 


dm/2 


-48 


14.87 


0.27 


bivalves 


dm/2 


-60 


14.82 


0.17 


bivalves 


dm/2 


-100 


14.88 


0.12 


bivalves 


dm/2 


-160 


14.90 


0.18 


bivalves 


dm/2 


-230 


14.91 


0.07 


bivalves 


dm/2 


-310 


15.09 


0.17 


bivalves 


dm/2 


-370 


15.25 


0.60 


bivalves 


dm/2 





12.3 


0.17 


bivalves 


ma 


-70 


12.6 


0.26 


bivalves 


ma 


-220 


12.6 


0.26 


bivalves 


ma 
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TABLE II. Original counting and geological data. 



Geologic time [M.Y.] 


Counting 


Adopted S.D. 


Organisms 


Data type 


Reference 


-300 


385 


4.3 


corals 


dy 


1 


2 b 




-330 


398 


6.1 


corals 


dy 


1 


28 




-390 


397.6 


1.8 


corals 


dy 


1 




! 


-425 


400 


6.1 


corals 


dy 


1 


2S 




-465 


412 


6.1 


corals 


dy 


[ 


_o 


I 





360 


3.1 


corals 


dy 


1 
1 


)S 

in 


1 


-390 


401 


6.1 


corals 


dy 


[ 


25 


! 


-395 


410 


4.5 


corals 


dy 


[ 


V 


! 


-440 


421 


6.4 


corals 


dy 


[ 


3C 


13 


-390 


30.63 


0.11 


corals 


dm 


I 
- 


29 


! 


-335 


30.2 


0.15 


corals 


dm 


1 

- 


}7 


1 


-480 


30.7 


0.23 


corals 


dm 


1 


32 




-395 


31.5 


0.45 


corals 


dm 


1 

- 


BC 


1 


-440 


32.4 


0.64 


corals 


dm 


| 


~ 


1 


-390 


13 


0.29 


corals 


ma 


[ 






-395 


13 


0.2 


corals 


ma 


I 
1 






-440 


13 


0.29 


corals 


ma 


[ 


3C 







359.3 


1.6 


bivalves 


da 


I 
1 


32 




-70 


375 


2.6 


bivalves 


da 




32 




-220 


371.6 


3.6 


bivalves 


da 




32 




-290 


383 


2.3 


bivalves 


da 




>— 




-350 


398 


1.8 


bivalves 


da 




>— 




-390 


405.5 


6.0 


bivalves 


da 


1 


32 







29.22 


0.08 


bivalves 


dm 


1 
1 


V> 




-10 


29.52 


0.08 


bivalves 


dm 


1 


!> 




-22 


29.42 


0.15 


bivalves 


dm 




™ 




-51 


29.7 


0.09 


bivalves 


dm 


1 


75 




-70 


29.85 


0.10 


bivalves 


dm 




)■) 




-220 


29.66 


0.15 


bivalves 


md 


1 
- 


32 


! 


-290 


30.16 


0.09 


bivalves 


dm 


1 

- 


32 


1 


-350 


30.37 


0.11 


bivalves 


dm 


[ 


32 


! 


-390 


30.35 


0.17 


bivalves 


dm 


I 
- 


32 


! 





12.35 


0.07 


bivalves 


ma 


1 
- 


y: 


! 


-70 


12.64 


0.11 


bivalves 


ma 


[ 

- 


32 




-220 


12.56 


0.08 


bivalves 


ma 


1 
- 


32 


1 


-1 


29.5 


0.08 


bivalves 


dm 


I 
- 


33 


! 


-5 


29.66 


0.15 


bivalves 


dm 


1 

- 


33 


1 


-30 


29.52 


0.08 


bivalves 


dm 


1 


33 


! 


-45 


29.72 


0.28 


bivalves 


dm 


I 


33 




-60 


29.7 


0.16 


bivalves 


dm 




33 




-100 


29.72 


0.05 


bivalves 


dm 


1 


33 




-170 


29.84 


0.3 


bivalves 


dm 


[ 


33 




-220 


29.76 


0.3 


bivalves 


dm 


1 


33 




-325 


30.13 


0.13 


bivalves 


dm 


1 


33 




-400 


30.5 


0.4 


bivalves 


dm 


1 


33 




-70 


29.65 


0.18 


bivalves 


dm 


1 
1 


71 




-70 


12.49 


0.02 


bivalves 


ma 


1 
1 






-320 


30.22 


0.40 


cephalopods 


dm 


1 






-410 


29.84 


0.23 


cephalopods 


dm 


I 
1 


32 




-310 


30.11 


0.35 


cephalopods 


dm 


1 


7f 




-395 


407.75 


3.2 


brachiopods 


da 


| 


I 




-425 


419 


4.5 


brachiopods 


da 




3C 




-395 


31.38 


0.32 


brachiopods 


dm 


1 


3C 




-425 


31.5 


0.45 


brachiopods 


dm 


1 


3C 




-395 


13 


0.15 


brachiopods 


ma 


1 


3C 




-425 


13 


0.2 


brachiopods 


ma 


1 


3C 




-540 


424 


6.4 


estromatolites 


da 


1 


35 




-150 


30 


0.64 


estromatolites 


dm 


I 


32 





G 



-510 31.56 0.75 estromatolites dm [32| 

-510 33 0.64 estromatolites dm [32 

-510 13 0.29 estromatolites ma [32 



TABLE III. Significance of the \ 2 test for all adjustments. 



Data source 


number of adjusted 
parameters 


X 2 significance 


Original 


5 


below 10% 


Original 


8 


below 10% 


Lambeck 


5 


17% 


Lambeck 


8 


85% 



TABLE IV. Adjusted parameters of the curves fitting 
Lambeck data, together with 95% confidence limits. 

Parameter best fit value confidence limits 

lower upper 

An© -1.3x10"^ -2.5x10"^ 1.6 xlO -3 

An, 9.3xl0~ 3 -4.1xl0~ 2 5.9xl0~ 2 

AQ 1.9xl0" 3 -4.1xl0~ 3 7.9xl0 -3 

£[10- n yr _1 16-4 12 20 

i[10- n yr _1 32.0 29 35 

§3- 6.5 xlO" 3 3.8xl0~ 3 9.2xl0~ 3 

period [10 6 yr] 486 370 680 

ip 0.27 -0.53 0.83 
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